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profiles  P (r)  ~ (1-r  /r  ) , modes  with  angular  dependence  Y „ (d,<p)  are 

o o JoItV 

unstable  provided  4 > 0 and  h < 1/(y~1)»  where  y is  the  ratio  of  specific 
heats.  For  large  4,  the  characteristic  growth  time  of  the  perturbations  varies 
as  4-ls  and  the  amplification  increases  exponentially  as  a function  of  4.  The 
radial  eigenfunctions  are  proportional  to  r , and  the  compressibility  and  vor- 
ticity  are  both  non-zero. 


I.  Introduction 

Over  a quarter  century  has  passed  since  the  discovery  that  the 
equations  of  one-dimensional  ideal  fluid  motion  possess  a variety  of 
self-similar  solutions.  These  were  found  independently  and  more  or 
less  simultaneously  by  Sedov  (1953),  Staniukovich  (1949),  Taylor  (1950) 
and  others,  stimulated  by  interest  in  nuclear  and  astrophysical  explo- 
sions and  in  the  general  properties  of  gas  dynamic  systems.  The  useful- 
ness of  these  solutions  is  twofold:  they  correctly  describe  1-D  flows 
at  late  times  when  the  details  of  initiation  or  preparation  have  been 
"forgotten,"  and  they  are  analytic,  or  at  least  reduce  the  solution  to 
quadratures. 

A particularly  useful  and  interesting  type  of  self-similar  motion 
is  that  known  as  uniform  or  homogeneous.  Its  characteristic  feature  is 
a velocity  field  which  is  proportional  to  the  distance  from  the  center 
of  symmetry.  Most  applications  have  been  to  problems  with  spherical 
symmetry,  such  as  supernova  explosions  (Keller,  1956),  laser  implosions 
(Kidder,  1976)  and  self-gravitating  clouds  (Sedov,  1959) . Closely  re- 
lated to  the  latter  are  cosmological  models  in  the  nonrelativistic 
limit  (Weinberg,  1972),  which  are  distinctive  by  virtue  of  being  pres- 
sureless and  unbounded. 

Sedov  (1959)  distinguishes  three  types  of  uniform  self-similar 
motion  in  an  ideal  gas.  In  type  I,  the  radius  varies  between  0 and  ». 

In  type  II,  it  varies  between  0 and  a finite  maximum  value,  correspon- 
ding to  a turning  point.  Type  III,  to  which  we  restrict  ourselves  in 
the  present  work,  has  radius  varying  between  a minimum  (at  a turning 


Note:  Manuscript  submitted  August  14,  1978. 


point)  and  °°.  We  follow  Keller  (1956)  in  using  Lagrangian  coordinates 
to  derive  in  §2.1  a two-parameter  family  of  solutions  including,  among 


others,  isothermal  and  uniform-density  models  as  special  cases.  There 
seems  at  first  glance  no  reason  why  these  should  not  be  stable,  and 
indeed  in  the  literature  where  applications  are  made  (e.g.,  Zel'dovich 
and  Raizer,  1966)  the  possibility  has  apparently  not  been  considered 
before.  We  go  on  in  §2.2,  however,  to  argue  that  instability  should 
arise  whenever  a condition  is  satisfied,  equivalent  to  the  presence  of 
an  entropy  density  which  decreases  outward.  Physically  the  mechanism  is 
identical  with  that  responsible  for  convective  instability  in  static 
stratified  media  when  the  temperature  decreases  in  the  upward  direction. 
Analysis  of  the  linearized  fluid  equations  in  §3  using  the  techniques 
developed  by  Bernstein  and  Book  (1978)  and  Book  (1978)  confirms  the 
existence  of  the  instability  and  yields  both  the  space  and  time  depen- 
dence of  the  perturbations  in  closed  form. 

As  was  previously  noted  by  Bernstein  and  Book  (1978)  and  Book  and 
Bernstein  (1978) , the  usual  definition  of  stability  is  inadequate  when 
applied  to  nonsteady  states,  since  the  time  dependence  of  the  perturba- 
tions is  in  general  not  exponential.  It  is  appropriate  to  call  a mode 
stable  (unstable)  if  the  ratio  of  the  perturbation  amplitude  to  that  of 
the  basic  state  vanishes  (diverges)  as  t -*■  °°.  We  find  that  in  the 
present  case  the  ratio  is  in  general  finite,  but  can  be  made  arbitrar- 
ily large  in  "unstable"  systems  by  choice  of  sufficiently  large  mode 
number 

A brief  discussion  of  the  results  in  §4  concludes  the  paper. 
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1.  The  basic  state 

1.1  Derivation  of  the  equations 

We  start  with  the  equations  of  ideal  hydrodynamics  which  in 
Lagrangian  variables  take  the  form 

p + v • Vp  = 0,  (2.1a) 

pv  + Vp  = 0 (2.1b) 

and 

(PP  Y)  =0  (2.1c) 

Here  dots  denote  time  derivatives  and  y in  (2.1c)  is  the  ratio  of 
specific  heats. 

In  a spherically  symmetric  system,  (2.1a,b)  become 

(2.2a) 

(2.2b) 

For  a motion  of  the  type  known  as  homogeneous  (or  uniform)  self-similar 
flow  (Sedov  1959) , the  position  R at  time  t of  a fluid  element  whose 
position  at  t = 0 was  r is  required  to  satisfy 

R = r f (t) , (2.3) 

where  f(0)  = 1 and  f(0)  = 0.  The  continuity  equation  (2.2a)  then  yields 

p(r,t)  = p (r) f 3,  (2.4) 

o 

and  hence  from  the  adiabatic  law  (2.1c) 

p(r,t)  = pQ(r)f  = s(r)p^  f 3Y  (2.5) 


PR"2  ~ (R2u) 


= 0, 


and 


pu  + 


3 


with  the  entropy  function  s arbitrary.  If  we  choose  the  initial 
density  profile  in  the  form 


P (r)  = (5(l-r2/r  2)  * 

O o 


where  0 and  k are  constants,  then  it  follows  from  (2.1b)  that 


and  f must  satisfy 


P (r)  = p(l-r2/r  2)  * +1, 
o o 


f f3Y~2  _ 2( K +l)p  _ ^-2 

0 r_2 


(2.6) 


(2.7) 


(2.8) 


We  will  use  the  separation  constant  t,  the  initial  radius  r , and  the 

o 

peak  mass  density  (3  to  rescale  t,  r and  p,  respectively.  In  these 
reduced  variables  we  have 


,,  2.  H 

PQ  = (1-r  ) , 


= 


(1-r2)  * +1 


o 2 ( H +1) 


(2.9) 


(2.10) 


and 


f f3Y'2  - 1. 


(2.11) 


A quadrature  can  be  performed  on  (2.11),  with  the  result 


f = 2&n  f 


(2.12) 


if  y = 1,  and 


(1  - f~a) 


(2.13) 


otherwise, where  a = 3(y-l).  If  y = 5/3,  (2.13)  can  be  integrated 

2 *5 

directly  to  give  f(t)  = ± (1  + t ) . For  other  values  of  y the  solu- 
tion is  most  conveniently  found  by  numerical  means.  At  large  |t|  when 
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f 0,  the  motion  asymptotically  approaches  free  streaming.  As  a 
function  of  the  anisentropicity  parameter  k , the  solutions  include 
the  cases  of  uniform  density  and  quadratic  pressure,  H = 0,  and  uni- 
form entropy  density,  x = l/(y-l) . If  we  exclude  singular  density 
profiles,  h is  restricted  to  0 s.  h < «®  . The  other  parameter  is  y, 
which  must  lie  in  the  range  1 £ y < *>. 

2.2.  Physical  mechanism  for  instability 

At  time  t a small  volume  AV  of  fluid  initially  located  at  radius  r 

“3  2 K+-1 

contains  a mass  Am  = pQ(r)AVf  , subjected  to  a pressure  p * (1-r  ) / 

3Y 

[2U+l)f  ].  Consider  two  such  fluid  elements  initially  at  radii  r^  and 

r2  > r^,  whose  volumes  are  related  by 

AV2/AV1  = (P1/P2)1/Y  = Kl-r12)/(l-r22)]U +1)/Y  . (2.14) 

The  compressional  energy  associated  with  these  elements  is 

Ep  = <P1AV1  + P2AV2)/(y-l) . (2.15) 

Their  kinetic  energy  calculated  from  the  expansion  or  contraction  of 
the  sphere  is 

(Am^u^2  + Am2u22)  . (2.16) 


Because  the  state  is  nonsteady  (*f  ^ 0) , the  elements  are  subject 

•• 

to  an  effective  gravitational  acceleration  g = rf,  and  therefore  have 
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a gravitational  potential  energy 


E = 'iff  (Am  r 2 + Am  r 2)  . (2.17) 

g 11  Z Z 


Now  let  the  two  fluid  elements  interchange  positions. 

By  (2.14),  they  contract  or  expand  so  as  to  satisfy  local  pressure 
balance  after  the  interchange.  Furthermore,  the  work  done  in  com- 
pressing one  is  just  balanced  by  that  done  by  the  expansion  of  the 
other,  so  the  compressional  energy  E^  afterwards  is  equal  to  E^.  The 
net  change  in  energy  is  then 


6E  = E"  + E.*  + E"  - E - E,  - E 
p k g p k g 


2 2 2 2 
= *5  [Am^u^  + Ami!u1  “ Am^u^  “ A®2U2  ^ 


+ *jff  [Am1r22  + Am2r12  - Am.^2  “ ^2*2  ^ 


= >j(f2  + ff)(r22  - r 2)  (Amx  - Am2) . 


(2.18) 


The  first  and  second  factors  in  the  last  member  of  (2.18)  are  strictly 
positive.  The  third  factor,  on  the  other  hand,  is  proportional  to 


2 . K +1 


“-•.•'■•“-'.''•(tv)  ’ 

. 2 . K +1-  H Y _ 

■*[-&)’  } 


(2.19) 


(2.19) 


This  expression  is  negative  for  H < l/(y-l) . In  this  case,  therefore. 


the  interchange  reduces  the  total  system  energy.  We  thus  anticipate 
that  an  instability  will  set  in,  characterized  by  "overturning"  of  the 
profiles,  such  as  is  typically  seen  in  convective  or  thermal  instabili- 
ties of  static  media  (Landau  and  Lifshitz,  1959) . 

When  H < l/(y-l) , 6E>  0,  in  which  case  no  instability  should 
arise.  The  marginal  case  just  corresponds  to  isentropic — more  properly, 
homentropic — states.  By  (2.5),  the  entropy  function  s satisfies 


The  stable  (unstable)  case  corresponds  to  outward  increasing 
(decreasing)  s(r).  Evidently  the  physical  picture  here  is  analogous  to 
that  arising  in  connection  with  instabilities  driven  by  a temperature 
inversion  in  media  with  a stratified  density.  Destabilization  takes  place 
owing  to  the  buoyancy  experienced  by  fluid  elements  in  the  nonuniform 
inertial  gravity  field.  It  is  therefore  purely  a consequence  of  the 
nonsteady  character  of  the  basic  state. 

3.  Analysis  of  the  perturbed  equations 

We  follow  Bernstein  and  Book  (1978)  and  Book  (1978)  in  obtaining 
linearized  equations  for  the  development  of  a small  perturbation  about 
the  solutions  of  §2.  For  simplicity  we  consider  only  expanding  states 
(t  > 0) . The  perturbed  displacement  £ satisfies  the  linearized  form 
of  (2.1b), 


(3.1) 


PI + p 1 i - - Vi + • 7R  p. 


Substituting  for  from  (2.3),  the  perturbed  density  from 


pi  = - p 7r  • k> 


(3.2) 


and  the  perturbed  pressure  from 


j = lE  p = - XillE  > 

•1  3p  M1  2(H  +1) 


2. 1+x  (1-Y) 


P VR  *i' 


(3.3) 


we  obtain  (writing  7 = 7 ) 


a+2  ••  r(l-r  )_ 

£ 2(*  +1)  ' ~ 


(Y-l)  r V • ( - V(  • r 


(3.4) 


Letting  a = 7 • 5 and  u>  = 7 x £,  we  have,  on  taking  the  divergence 
and  curl  of  (3.4), 


*0+2  •• 
f a = 


2(h  +1) 


7 • [(1-r  )7o]  - yr  • 7a 


- (3y-2)o  + r • 7 X a) 


and 


,a+2 


to  = [ 1 - -^rr|7o  X r + to 
~ I x +1 1 ~ 


(3.5) 


(3.6? 


We  look  for  solutions  of  (3.5)-(3.6),  assuming  £ is  separable  into  a 
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product  of  a function  of  position  and  a factor  T(t)  satisfying 


f 


ot+2 


T = pT, 


(3.7) 


P constant.  We  further  assume  separation  of  the  angular  and  radial 
dependence  by  writing 


o (r)  = o (r)  Y (6,40  . (3.8) 

In  (3.5)  to  appears  only  in  the  form  r • V x to,  for  which  an  expres- 

^ rv 

sion  in  terms  of  o can  be  derived  from  (3.6)  and  (3.8): 

(p-l)r  • V x w = K [2r*Vo+rr: VVo-r2V2o] 

= a [1- (y-1)  *]  «.(«.+  1)/(h+1).  (3.9) 


Substitution  in  (3.5)  yields  a second-order  equation  for  the  radial 
factor  a (r) , 


Y f~  1-r2  _d_ 

!(*+l)  2 dr 

L r 


, 2 da.  1-r 
(r  2 


ZU+l)o  - 2r 


do 

dr 


yr 


do  ( K +1-  H y ) l 

dr  + |_  (p-1)  ( h 


& ( &+1) 
+1) 


i o = 


- 3y  - U+2  0 = 0.  (3.10) 


Rewriting  this  equation  by  means  of  the  substitutions  o = r y 


and  x = r , we  obtain  the  hypergeometric  equation 
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where 


x(l-x)y"  + [c  - (a+b+l)x]  yy  - aby  = 0, 


(3.11) 


||  = *|*+A  + 5/2±  [(*  + l + 5/2) 2 - 4K]1^, 


c = l + 3/2. 


(3.12a,b) 


(3.12c) 


( H +2)1,  H +1 

2 2'y 


-2  + 3y 


&U+1) 

(M-l 


■imx  +D-w  Yn 

-1)  ( x +1)  J.  (3*13) 


The  solution  of  (3.11)  which  is  finite  at  the  origin  is  the  hyper- 
geometric function  y = 2F^(a,b;c;x) . 

The  boundary  condition  is  found  from  the  requirement  that  the 
perturbed  pressure  vanish  on  the  interface  with  the  vacuum.  Since  the 
unperturbed  pressure  already  vanishes  there,  it  follows  from  (3.3)  that 
y need  only  be  finite  at  x = 1.  The  analytic  continuation  formulas 
(e.g.,  Abramowitz  and  Stegun  1968)  contain  a term  that  diverges  as 
(1-x)  ^ K unless  a or  b is  a nonpositive  integer.  Thus  we  must  have 


n = hi  H +1+5/2  - tOt+fl+5/2)2  - 4K]*5}  , 


(3.14) 


n=0,  1,  2,  •••  . Treating  n and  K as  analytic  functions  of  y and 
differentiating  (3.14)  yields  9n/3y  < 0.  Hence  the  fastest  growth 
(Largest  y>l)  corresponds  to  the  smallest  value  of  n,  viz.,  n = 0, 
which  implies  K = 0.  Solving  for  y,  we  finally  obtain  the  dispersion 


relation 


L 


U 1 = [ +2)+(*  +1)  (3y-1)  1 

M 2 ( H +1) 


± Uy&(*  +2)+(*  +1)  (3y-l)  ]2+4Jm+l)  (H  +1-*  y)  }** 

2(  K +1) 


(3.15) 


For  the  upper  branch,  u>l  for  all  £>0,  provided  H<l/(y-l).  The  latter 
is  precisely  the  condition  derived  from  the  energetic  argument  of 
S2.2. 

Returning  to  (3.7),  we  find  that,  provided  y>l,  the  time  depend- 
ence can  likewise  be  expressed  in  terms  of  hypergeometric  functions  in 
the  form  (Bernstein  and  Book,  1978) . 


T(t)  * T(0)  3-  (t)  + f(0)  S (t). 


(3.16) 


Here, 


■ 2P11I  ♦ TT-  i * Hr  ? I-'-0).  13.17,) 


i It)  - [£  U-£-“))  Vi  Ij  l-f'“).  (3.17b) 


r 2 , J. 

and  A = l(o+2)  - 8yaJ  . For  late  times  (large  f) , the  analytic 


continuation  formulas  yield 


3-<t) 

Sit) 


r(l/2)r(l/g)  f 


T [ 1/4+ ( 2+A) /4a ] T [ 1/4+ ( 2-A) /4a ] ' 


(3.18a) 


(2/a) " r(3/2)r(l/g)  f 


r[3/4+(2+A)/4o]r[  3/4+(2-A) /4a] 


(3.18b) 
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The  numerical  coefficients  in  (3.18)  grow  exponentially  with  y for  y >>  1. 
The  case  of  y=l  is  very  similar,  except  that  confluent  hypergeometric 
functions  replace  2Fi'  as  observed  by  Bernstein  and  Book  (1978)  and  Boo!: 

and  Bernstein  (1978),  and  (3.18)  is  replaced  by  expressions  proportional 

„ ,,  ,,  (P-U/2 

to  f (to*  f ) 

4.  Discussion 

We  have  seen  on  energetic  grounds  that  a certain  class  of  spheri- 
cal ideal  gas  expansions  can  be  expected  to  be  unstable  whenever  the 
gradient  of  the  entropy  density  decreases  with  increasing  r.  Detailed 
analysis  of  the  linear  perturbations  about  these  nonsteady  basic  states 
confirms  this  prediction,  provided  we  appropriately  generalize  the 
usual  definition  of  instability.  Somewhat  surprisingly,  the  solutions 
fall  out  exactly  without  recourse  to  numerical  approximations,  owing  to 
the  separability  of  the  linearized  equations. 

As  noted  in  §1,  what  matters  in  determining  the  stability  of  a 
time-dependent  motion  is  the  relative  size  of  the  perturbations.  By 
(3.18),  the  latter  vary  asymptotically  like  the  unperturbed  radius. 

At  early  times,  however,  when  f differs  substantially  from  zero,  the 
perturbations  can  be  amplified  dramatically.  If  y >>  1,  they  grow 
approximately  exponentially  for  t ^ 1,  experiencing  ~ y e-foldings. 

The  total  amplification  and  the  time  required  to  approach  the  asymp- 
totic state  in  which  they  "freeze  out"  both  increase  with  y.  As  y -*•  1, 
both  the  total  amplification  and  the  time  required  to  approach  satura- 
tion diverge  (Bernstein  and  Book,  1978) . Since  y increases  with 
increasing  l,  decreasing  x , and  decreasing  y,  ell  of  these  trends  tend 
to  enhance  instability. 


12 


Note  that  as  £ •+  00 , y diverges.  This  implies  that  the  problem  is 
not  mathematically  well-posed.  In  any  real  physical  system,  dissipa- 
tive phenomena  related  to  viscosity,  thermal  conduction,  radiation, 
etc.,  set  an  upper  limit  on  the  mode  number  for  which  the  ideal  fluid 
model  is  valid.  For  shorter-wave length  disturbances  than  this,  not 
only  the  detailed  perturbation  analysis,  but  the  whole  physical  picture 
must  be  drastically  different. 

The  perturbations  studied  here  have  radial  dependence  which  peaks 

at  r = rQ.  They  therefore  should  be  most  readily  observable  as  an 

enhanced  mixing  or  turbulence  near  the  periphery  of  the  expanding 

cloud.  Since  the  instability  is  controlled  by  the  sign  of  the  entropy 

gradient,  it  seems  likely  that  the  nonlinear  limit  to  which  it  tends 

is  characterized  by  ds/dr  ^ 0,  0 £ r £ r . Whether  this  limit  is 

* o 

actually  attained  is  beyond  the  scope  of  the  present  work. 

Another,  perhaps  more  important,  question  remains  unanswered. 
Uniform  self-similar  motion  is  an  analytically  convenient  model  used 
to  approximate  real  flows.  To  what  extent  is  the  instability  treated 
here  associated  with  the  latter,  to  what  extent  an  artifact  of  the 
model?  The  present  paper  can  of  course  provide  no  rigorous  answer. 
Nonetheless,  it  seems  physically  plausible  that  for  flows  sufficiently 
close  to  uniform  expansion,  the  results  of  the  present  analysis 

must  be  applicable.  Even  for  nonuniform  motions,  either  analytically 
or  numerically  described,  the  energetic  argument  of  §2.2  can  be  em- 
ployed and  should  again  correctly  predict  the  presence  or  absence  of 
instability. 
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